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Abstract 

We consider the fractional generalizations of equation that defines the medium 
mass. We prove that the fractional integrals can be used to describe the media with 
noninteger mass dimensions. Using fractional integrals, we derive the fractional gener- 
alization of the Chapman-Kolmogorov equation (Smolukhovski equation). In this pa- 
per fractional Fokker-Planck equation for fractal media is derived from the fractional 
Chapman-Kolmogorov equation. Using the Fourier transform, we get the Fokker- 
Planck-Zaslavsky equations that have fractional coordinate derivatives. The Fokker- 
Planck equation for the fractal media is an equation with fractional derivatives in the 
dual space. 

The application of fractals in physics is far ranging, from the dynamics of 
fluids in porous media to resistivity networks in electronics. The cornerstone of 
fractals is the meaning of dimension, specifically the fractal dimension. Fractal 
dimension can be best calculated by box counting methods which means drawing 
a box of size R and counting the mass inside. Fractal models of media are 
enjoying considerable popularity. This is due in part to the relatively small 
number of parameters that can define a fractal medium of great complexity and 
rich structure. Derivatives and integrals of fractional order have found many 
applications in recent studies of scaling phenomena. It is interesting to use 
fractional integration to consider the fractional generalizations of the equations 
that describe fractal media. We use the fractional integrals to describe fractal 
media with noninteger mass dimensions. The fractional integrals can be used not 
only to calculate the mass dimensions of fractal media. Fractional integration 
can be used to describe the dynamical processes in the fractal media. Using 
fractional integrals, we derive the fractional generalization of the Chapman- 
Kolmogorov equation (Smolukhovski equation). The fractional Fokker-Planck 
equation for fractal media is derived from the fractional Chapman-Kolmogorov 
equation. Using the Fourier transform, we get the Fokker-Planck-Zaslavsky 
equations that have coordinate fractional derivatives. 



1 Introduction 



The application of fractals in physics |H 12 El H] is far ranging, from the dynamics of fluids 
in porous media to resistivity networks in electronics. Derivatives and integrals of fractional 
order have found many applications in recent studies of scaling phenomena El El El 13 El E] ■ 

The natural questions arise: What could be the physical meaning of the fractional inte- 
gration? This physical meaning can be following: the fractional integration can be considered 
as an integration in some fractional space. In order to use this interpretation we must define 
a fractional space. The first interpretation of the fractional space is connected with frac- 
tional dimension space. The fractional dimension interpretation follows from the formulas 
for dimensional regularizations. If we use the well-known formulas for dimensional regular- 
izations PU]> then we get that the fractional integration can be considered as an integration 
in the fractional dimension space ^1] up to the numerical factor r(a/2)/[27r a//2 r(a)]. This 
interpretation was suggested in Ref. jTT]. 

In this paper we use the second physical interpretation of the coordinate fractional inte- 
gration. This interpretation follows from the fractional measure of space ^T] that is used in 
the fractional integrals. We consider the mass fractional dimension and the fractional gener- 
alizations of the equation that defines the medium mass. We use the fractional integrals to 
describe fractal media with noninteger mass dimensions [12]. Using fractional integrals, we 
derive the fractional generalization of Chapman-Kolmogorov equation (Smolukhovski equa- 
tion). In this paper the fractional generalization of Fokker-Planck equation is derived from 
the fractional Chapman-Kolmogorov equation. Using the Fourier transform, we derive the 
Fokker-Planck-Zaslavsky equations El that have fractional coordinate derivatives. The 
Fokker-Planck equation for the fractal media is an equation with fractional derivatives in 
the dual space. 

In Sec. 2, the definition of the mass dimension and the fractional generalization of the 
media mass equation are considered. In Sec. 3, we derive the fractional generalization of 
the average values equation. In Sec. 4, the fractional Chapman-Kolmogorov equation is 
derived by using fractional integration. In Sec. 5, the fractional Fokker-Planck equation for 
the fractal media is derived from the suggested fractional Chapman-Kolmogorov equation. 
In Sec. 5, we use only the fractional power coordinate derivatives. The fractional derivatives 
are used only in Sec. 7. In Sec. 6, the stationary solution of the Fokker-Planck equation for 
fractal media are derived. In Sec. 7, we consider Fokker-Planck-Zaslavsky equations with 
fractional coordinate derivatives. We derived these equations by the Fourier transform of 
the Fokker-Planck equations for fractal media. 

2 Mass Fractal Media and Fractional Integral 

Equations that define the fractal dimensions have the passage to the limit. This passage 
makes difficult the practical application to the real fractal media. The other dimensions, 
which can be easy calculated from the experimental data, are used in empirical investigations. 
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For example, the mass fractal dimension ^Sl Ej can be easy calculated. The properties of 
the fractal media like mass obeys a power law relation 



M(R) = kR Um , (1) 

where M is the mass of fractal medium, R is a box size (or a sphere radius), and D m is a 
mass fractal dimension. The amount of mass of a medium inside a box of size R has a power 
law relation ([[J. 

Fractal dimension can be best calculated by box counting which means drawing a box 
of size R and counting the mass inside. To calculate the mass fractal dimension, take the 
logarithm of both sides of Eq. (JTJ): D m = ln(M)/ln(R). Log-log plot of M and R gives us 
the slope D m , the fractal dimension. When we graph ln(M) as a function of ln(R) we get a 
value of about D m which is the fractal dimension of fractal media. 

Let us prove that power law relation (0) can be naturally derived by using the fractional 
integral approach [TI]. Let us derive the fractional generalization of the equation, which 
defines the mass of medium that is distributed on the line —(R/2) < x < (R/2) and has the 
form 

rR/2 

Mi(R) = / p{x)dx, (2) 



-R/2 

where p(x) G This equation can be written in the following form: 

rR/2 



-R/2 

Using the notations ^3] for fractional integrals 



Mi(R)= / p(x)dx+ / p(x)dx. (3) 

J -R/2 Jy 



fra w \ 1 [ V p(x)dx 



xV- a 



MM 1 ,m P[x)dx 



where —(R/2) < y < (R/2), we can rewrite Eq. (jSJ) in the form 

M 1 (R) = (I 1 R/2+ p)(y) + (I 1 R/2 _p)(y). 
The fractional generalization of this equation is defined by 

M a (R) = (I a R/2+P )(y) + (I^p)(y), (4) 

where p(x) G L p (R 1 ), and 1 < p < 1/a. The initial points in Eq. (J3J) are set to zero (y = 0). 
As the result, we can rewrite Eq. (j3J) in an equivalent form 

-R/2 



M a (R) = — / ip(x) + p(-x)]x a - l dx. 
1 a jo 
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In order to describe the fractal media, we suggest to use the fractional continuous medium 
model ^2]- We proposes ^2] to describe fractal media by a "fractional" continuous medium 
model where all characteristics are defined everywhere in the volume but they follow some 
generalized equations which are derived by using fractional integrals. 

If we use Eq. (jlj) and use the constant density distribution (p(x) = p = const) for 
homogeneous fractal media, then we get 



We have that the mass fractal dimension D m = a, and k = po2 1 ~ a /[aT(a)] in Eq. The 
fractional generalization of the equations for fractal media in the three-dimensional space 
are considered in Ref. ^2] and in the Appendix of this paper. 

The fractional integrals can be used to describe fractal media with noninteger mass di- 
mensions j!2j . The fractional integrals can be used not only to calculate the mass dimensions 
of fractal media. Fractional integration can be used to describe the dynamical processes in 
the fractal media. Using the fractional integrals ^Tj, we can derive the fractional general- 
izations of the Chapman-Kolmogorov and Fokker-Planck equations. 

3 Fractional Average Values 

In order to derive the fractional analog of average value we consider the fractional integral 
for the function A(x). The usual average value 



Po2 



l-Q 



M a (R) 




R 



i a: 




(5) 



can be written in the form 





Using the following notations for the fractional integrals, 




(7) 





the average value © can be rewritten in the form 



<A>, (y) = (I 1 + Ap)(y) + (I 1 _Ap)(y). 
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The fractional generalization of this equation can be defined by the equation 

< A > a (y) = (ItAp)(y) + {I a _Ap){y). (9) 
Note that Eq. Q can be written [TT] in the equivalent form, 

f+OO 

<A> a (y)= [(Ap)(y-x) + (Ap)(y + x)}dfi a (x), (10) 
Jo 

where we use 

|x| a_1 <ix dx a 

dpi a (x) - 



r(a) aT{a) 

Here we use the following notation for fractional power of coordinates: 

x a = (3(x)(x) a = sgn(x)\x\ a } (11) 

where /3(x) = [sgn^)]^ 1 . The function sgn(x) is equal to +1 for x > 0, and this function 
is equal to —1 for x < 0. 

In order to have the symmetric limits of the integral we consider average value (fTUj) in 
the form 

| f+OO 

<A> a (y) = - [(Ap)(y-x) + (Ap)(y + x)]dfi a (x). (12) 



2 

If a = 1, then we have the usual equation for the average value. 

Let us introduce some notations in order to simplify the equation for the fractional 
average value. We use the integral operators I" that is defined by 



1 



-oc 



!*/(*) = o/ f(x) + f(-x) dfi a (x). (13) 

/ J—oo L J 

In this case, fractional generalization of the average value ()12j) can be rewritten in the form 

< A> a =r x A(x)p(x). (14) 

We will use the initial points in Eq. (J3J) are set to zero (y = 0). Note that the fractional 
normalization condition ^T] is a special case of this definition of average values < 1 > a = 1. 



4 Fractional Chap man- Kolmogorov Equation 

The Fokker-Planck equation can be derived from so-called Chapman-Kolmogorov equation 
(also known as Smolukhovski equation) [TBI El CHI • This equation may be interpreted as 
the condition of consistency of distribution functions of different orders. 

Kolmogorov [THl [T7| derived a kinetic equation using a special scheme and conditions that 
are important for kinetics. Let W(x, t; Xq, to) be a probability density of having a particle at 
the position x at time t if the particle was at the position xq at time t Q < t. 
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Denote by p(x, t) the distribution functions for the given time t. Let us consider two 
well-known identities 

/+oo r+oo 
dx' W(x,t\x',t')p(x',t'), / p(x,t) = l. (15) 
-oo J — oo 

Using notation (fTT?|) . we can rewrite these equations in the form 

p(x, t) = % W(x, t\x', t')p{x', t'), ll p(x, t) = 1. 

Therefore we can get the fractional generalization of Eqs. (fTK|) in the form 

p(x, t) = l a x , W{x, t\x\ t')p(x', t'). (16) 

This equation is the definition of conditional distribution function W(x, t\x', t') for the fractal 
medium referring to different time instants. The normalization conditions for the functions 
W(x, t\x', if) and p(x,t) have the form 

I«W(x,t\x',t') = l, l a xP (x,t) = l. (17) 

Substituting into the right-hand side of Eq. (jlfi)) the value of p(x', t') expressed via the 
distribution p(x ,t ) at an earlier time, 

p{x', t') = 1° W{x', t'\x , t )p{x , t ), (18) 

we obtain the integral relation which includes the intermediate point x', 

p(x, t) = % l a xo W(x, t\x', t')W(x', t'\x , to)p(x , t ). (19) 

Using Eq. (|T^|) . and Eq. ()16|) in the form 

p(x, t) = W(x, t\x , t )p(x , t ), (20) 

we can derive a closed equation for transition probabilities 

l a X0 W(x,t\xo,to)p(xo,to) = i a x , l a X0 W(x,t\x',t')W(x',t'\xo,to)p(xo,to). 

Since the equation holds for arbitrary p(xo,t ), we may equate the integrand. As the result, 
we get the fractional Chapman- Kolmogorov equation 

W(x, t\x , to) = I x , W(x, t\x', t')W(x', t'\x , t ). (21) 

This equation describes the Markov-type process in the fractal medium that is described 
by the continuous medium model [TS]. Using this fractional equation we can derive the 
fractional Fokker-Plank equation. 



5 Fractional Fokker-Planck Equation from 
Fractional Chap man- Kolmogorov Equation 

Let us consider the fractional generalization average value (|14j) for the function A(x), 

< A> a =i% A(x)p(x,t). (22) 
We will use the initial points are set to zero (y — 0). Using Eq. (JT5J1 in the form 

p{x, t) = l a XQ W(x, t\x , t )p(x , t ), (23) 

we get the average value 

< A > a = I£ A{x) l a X0 W(x,t\xo,t )p(x ,t ). 
We can rewrite this equation in the form 

< A > a = i« p(x , t ) I a x A(x) W(x, t\x , t ). (24) 

Let us consider the Taylor expansion of function A(x). If we use this expansion in the 
form 

A( X ) = A( X0 + Ax) = A M + A* + JA.C + .... (25) 

where Ax = x — Xq, then the fractional integration by parts of Eq. (J24j) cannot be realized 
in the simple form. It can be assumed jTTj that the function A has the form A = A(x a ). 
The Taylor expansion can be used in the form 

A(a-) = A(x% + Ax a ) = A«) + (^l) x Ax^ + l{^^)j^ a ) 2 + .~, (26) 

where Ax a and x a = sgn(x)\x\ a is defined by Eq. (II Here and later, we use 

the following notations for derivatives of fractional power of coordinates: 

d \x\^ d (2?) 



dx a a dx 

Note that these derivatives are the derivatives of fractal order in the dual space ("momentum 
representation"). The expansion ()2(j)l is equivalent to the usual Taylor expansion. In order 
to prove this proposition, we consider the Taylor expansion for x a that has the form 

Ax a = x a - x« = Ax + 1(^) (Ax) 2 + ... . (28) 

v ox J x o 2 y (oxy j x " 



Substituting Eq. (J2~H|) into Eq. (j2T?j) . we get the usual Taylor expansion (|23|) . If we use the 
usual Taylor expansion (|2*3j) . then the integration by parts in Eq. (J2~4"j) is more complicated. 
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If we use the Taylor expansion for the fractional power of coordinates, then the integration 
by parts in Eq. ()24j) can be realized in the simple form, 



dx a J-oo aiT (a) dx a 

+00 r+00 dx a a ,dB(x) 



= (B(x)A(x)) - / ——A (x' , 
v /-00 j_ OQ a r(a) ox 

Substituting Eq. ((23) in Eq. (1241 . we get 



< A > a = I a XQ A{x%)p(x ,t )i° W(x,t\x ,t )+ 

+ \t (^^)^°'« (A^) a WM*o,*o) + - (29) 
Let us introduce the following functions: 

P n (x , t, t ) = I" (Ax a ) n W(x, t\x , t ). (30) 
Using (f3T)|) and fTTj) . we can rewrite Eq. ()29j) in the form 



< A > a = 1% A(x%)p(x Q ,t Q ) +i a xo (-^^) xo p(x ,t )P l (x ,t J t ) + 
2 x ° V (dx a ) 



+ kc (^?) p(x ,t )P 2 (x ,t,t ) + ... • (31) 

XQ 



Using (j22j) in the form 

<A> a =l a Xo A(x%)p(x ,t), 
for the left hand side of Eq. (J3~T|) . we can rewrite this equation 

to ^(^o)(p(a;o,t) -p(x ,t ; 



= t (^?),/(^,to)P 1 (xo,t,t ) + it {j^l) xo PMP 2 (xo,t,t ) + ... . (32) 

The final step is an assumption that is usually called Kolmogorov condition ^SJ E] ■ We 
may assume that the following finite limits exist: 

,. Pi(x,t,t ) 

hm = a(x, tn), 

At^o At v 1 

V P 2(x,t,t ) 

ii£b At =6(Mo) ' 
8 



and 

P n (xt,t ) = 

At^O At 

where n > 3. It is due to the Kolmogorov conditions that irreversibility appears at the final 
equation. Multiplying both sides of Eq. (J3~2"|) by 1/At and consider the limit, we obtain 



dp(x , t) 



dt /to 



Integrating by parts, we obtain 

Z^**M*,t) = -t:A<.*) e «'M* t » , (33) 

igW^ tWj|t) = i . A(0 W(M» (34) 



Here we use 



lim p(x, t) = 0. 

x— >±oo 



Then we have the equation 

fa Af^f dpfat) , d(p(x,t)a(x)) _ l d 2 (p(x,t)b(x)) ^ _ 

Since the function A = A(x a ) is an arbitrary function, then we have the following equation 



dp(x,t) d(p(x,t)a(x,t)) l d 2 (p(x,t)b(x,t)) _ 



<9t <9a; a 2 (<9r 

that is the Fokker-Planck equation for fractal media. This equation is derived from the 
fractional generalization of the average value and fractional normalization condition, which 
use the fractional integrals. The Fourier transform of this equation has the fractional deriva- 
tives. We consider the equations with fractional derivative in Section 7. Therefore this 
equation can be called the fractional Fokker-Planck equation. The Fokker-Planck equation 
with fractional derivatives we call Fokker-Planck-Zaslavsky equation, since this equation was 
suggested by Zaslavsky 121 El- 
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6 Stationary Solutions 

In this section, we consider simplest solutions of the Fokker-Planck equations for the fractal 
media. Let us consider the stationary Fokker-Planck equation for fractal media 

d(p(x,t)a(x,t)) _ l d 2 (p(x,t)b(x,t)) = 

dx a 2 (dx a ) 2 ' 1 ' 

This equation can be rewritten in the form 

^W)- ^'" )^ (37) 
Obviously, we get the relation 

p(x , t)a{x , t) -\ ?^mEM = cmlst . (38) 



Supposing that the constant is equal to zero, we get 

d(p(x,t)b(x,t)) 2a(x,t) 



dx a b(x, t) 

or, in an equivalent form 



{p{x,t)b{x,t)), (39) 



(40) 



d ln(p(x,t)b(x,t)) _ 2a(x,t) 
dx a ~~ b(x,t) 

The solution of this equation is 

ln(p(x, t)b(x, t)) = f 2 w ! dx a + const. (41) 

As the result, we have the following solution of the stationary Fokker-Planck equation for 
fractal media 

B^r** 2 /SiK (42) 

where the coefficient iV is defined by the normalization condition. 

Let us consider the special cases of this solution. 
(1) If a(x) = k and b(x) = —D, then Fokker-Planck equation for fractal media has the form 

dp(x,t) u dp(x,t) D d 2 p(x,t) _ n 
The stationary solution of this equation is the following distribution function: 



a 



I 1h\x\ 

p(x,t) = N 1 exp[ ^-J. (44) 



D 
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(2) If a(x) = k\x\P and b(x) = —D, then Fokker-Planck equation for fractal media has the 
form 

dp(x,t) u d\x\^p(x,t) D d 2 p(x,t) _ n 

~^r~ + k dx~« + 2"7^F ~ [ } 

The stationary solution of this equation is the following distribution function: 

p(x, t) = N 2 exp (-^—LL^ J ■ (46) 

If P = a, we get 

p(x,t) = N 2 exp (-^\x\ 2( *j . (47) 

If a + (5 = 2, we have 

p(x, t) = N 2 exp ^-^^ 2 J • (48) 

(3) If a(x) = k\xf and b(x) = — -D|x| -7 , then the stationary solution of the Fokker-Planck 
equation for fractal media is the following distribution function: 

p(x,t) = N 3 \x\iexp (- {a+ 2 pl i)D \x\ a+P+ ^ ■ (49) 

(4) If the coefficient a(x) is defined by 

, , dUix) \x\ l - a dU(x) 

aix) = — = — , 

ox a a ox 

and b = —D, then we get the following stationary solution: 

U(x)\ 



p(x, t) = N^exp 



D 



Let us consider the Fokker-Planck equation for fractal media ()35|) with a(x) = k\x\ a and 
b = —D. The stationary solution (|4T|) has the form 



p(x) 



k \ 1/2 _,.,2, 
7rD 



kx 2a /D 



The general solution can be represented in the form 



+°° / b i 

pOm) = E \l^^e- kx2a/DH n(x«V^ )e -^A n , 



n=0 

where 



An = V 2^! ^ P( x >°) H n(xyk/D)- 

Note that 

< x a (t)x a (0) > a = ^e~ kt . 
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7 Fokker-Planck-Zaslavsky Equations 



The fractional generalization of the Fokker-Planck equations was suggested by Zaslavsky 
in Ref. j^j. These equations have fractional derivatives and can be used to describe the 
dynamical processes on fractals [21 EH E] • The fractals can be realized in nature as fractal 
media. Therefore it is interesting to derive Fokker-Planck-Zaslavsky equation as an equation 
for fractal media [T2*] . 

We derive the Fokker-Planck-Zaslavsky equations that have coordinate fractional deriva- 
tives in the dual space. The Fokker-Planck equations for fractal media are equations with 
fractional derivatives in the dual space ("momentum representation"). Let us use Eq. (|35|) 
and the Fourier transform T of the coordinates that leads us to the Reisz derivatives. 

Using Eq. (|27|). we can rewrite the Fokker-Planck equation for fractal media in the 
form with usual coordinate derivatives d x = d/dx: 

| r |l— a (1 _ „,\| T |l-2a I ™ 1 2— 2a: 

dtpit, x) + J±!_a,(o(x)p(t, x)) - { - ^ ^.(6(x)p(t, x)) - m- 5 -flg(6(x)p(t, x)) = 0. 

a 2a z 2a z 

Let us consider two special cases, (1) the coefficients a(x) and b(x) are constant, (2) the 
coefficients a(x) and b(x) are connected by the equation a(x) = (l/2)db(x)/dx a . 

(1) If the functions a(x) and b(x) do not have the coordinate dependence (a = k and 
b = —D), then we get the equation 

dp(x,t) u dp(x,t) D d 2 p(x,t) _ n 

~^ + k ^x^ + 2"7^F" U - (5Uj 

This equation can be rewritten in the form 

o / x k\x\ l ~ a n , x D(l — a)|x| 1_2a , . D\x\ 2 - 2a ^ , s 

d t p{t, x) + d x p(t, x) — d x p{t, x) — — d x p[t, x) = 0. 

a 2a z 2a z 

The Fourier transform T of the coordinates \x\^ leads us to the Reisz derivatives ^H] 
p _ 1 ( 1 d [v f(z)dz -1 d [+™ f(z)dz 



y 2cos{<Kf3/2)^T{l-f3)dyJ-oo{y-zy T(l - (3) dy J y (z - yf 
by the equations 

F[\xfp{t,x)] = -nPp{t,y), 

where p(t, y) = T[p{t, x)\. As the result, we have the fractional Fokker-Planck equation with 
fractional derivatives in the dual space 

d t p(t,y) + ^Dj-fy/Kf,!/)) + D( 3~ a) Di- 2a (»yp(t,y)) + £- 2 V 2 y - 2a (y 2 p(t, y)) = (51) 
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for the function p(t,y). Here we use jF[<9"p(t, x)] = (—iy)p(t,y). Equation (JoTj) is a dual 
form of Fokker-Planck equation (|5()jl . As the result, we have fractional Fokker-Planck in the 
dual space. 



(2) After using the relation 



a{x) = (52) 



we get the diffusion equation (J33|) in the form 

with diffusion coefficient _D(x) = b(x). Here derivatives d/dx a are defined by Eq. (|27|). 
Using Eq. (}2Tj) . we have 



dp(x,t) \x 



l-a 



t -d x (\x\ 1 - a D(x)d xP (x,t)) =0. 



^ 2a 2 

If the diffusion coefficient D(x) is described by the power law 

\x\ 1 - a D{x) = 7 |x| /3 , 

then we have the diffusion equation in the form 

d tP (x,t) - ^L\x\ l -<*d x (\x\Pd x p(x,t)) = 0. 

The Fourier transform T of the coordinates leads us the following fractional Fokker-Planck 
equation with the Reisz fractional derivatives 



d t p(y,t) - ^Dj- Q (yDjyp(y,t)) = 0. 



If /3 = 0, then we have the fractional kinetic equation 

d t p(y,t)=r>l- a (A(y)p(y,t)), (53) 

where A(y) = (7 /2a 2 )y 2 . As the result, we get the equations with fractional derivatives that 
have Fokker-Planck-Zaslavsky form j2j |3j in the dual space. 
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8 Conclusion 



The application of fractals in physics is far ranging, from the dynamics of fluids in porous 
media to resistivity networks in electronics. The cornerstone of fractals is the meaning of 
dimension, specifically the fractal dimension. Fractal dimension can be best calculated by 
box counting method which means drawing a box of size R and counting the mass inside. 
When we graph ln(M) as a function of hx(R) we get a value of about D m which is the 
fractal dimension of fractal medium. The mass fractal dimensions of the media can be easy 
measured by experiments. The experimental determination of the mass fractal dimension 
can be realized by the usual box counting methods. The mass dimensions were measured 
for porous materials [1211201, polymers [21], and colloid agregates [22J. Fractal models of 
porous media are enjoying considerable popularity [231 1211 123 123 121] • This is due in part to 
the relatively small number of parameters that can define a fractal porous medium of great 
complexity and rich structure. 

We suppose that the concept of fractional integration provides an alternative approach 
to describe the fractal media. In this paper we prove that fractal media can be described by 
using the fractional integrals. We consider the fractional generalizations of the equation that 
defines the medium mass. This fractional generalization allows us to realize the consistent 
description of the fractal mass dimension of the media J21- The fractional integrals can be 
used in order to describe the dynamical processes in the fractal media. Fractional integration 
approach is potentially more useful for physics of fractal media than traditional methods that 
use the integer integration. In this paper fractional generalizations of the average value is 
derived. Using fractional integrals, we derive the fractional generalization of the Chapman- 
Kolmogorov equation. In this paper we derive the Fokker-Planck equations for the fractal 
media from the suggested fractional Chapman-Kolmogorov equation. Using the Fourier 
transform, we get the Fokker-Planck-Zaslavsky equations that have coordinate fractional 
derivatives in the dual space ("momentum representation"). 

Note that the Liouville equation is a cornerstone of the statistical mechanics [23 123 133 ■ 
The fractional generalization of the Liouville equation was suggested in Ref. [TI]. The frac- 
tional generalization of Liouville equation allows us to derive the fractional generalization of 
the Bogoliubov equations, Fokker-Planck equation, Enskog transport equation, and hydro- 
dynamic equations. Using fractional analog of the Liouville equation ^T] , we can derive the 
fractional Fokker-Planck equation for the phase space. 

Using the methods suggested in Refs. [33 133 133] , we can realize the Weyl quantization 
of the suggested fractional equations. 
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Appendix 



In many problems the real fractal structure of matter can be disregarded and the medium can 
be replaced by some "fractional" continuous mathematical model In order to describe 
the media with noninteger mass dimension, we must use the fractional calculus. Smoothing 
of the microscopic characteristics over the physically infinitesimal volume, we transform the 
initial fractal medium into "fractional" continuous model that uses the fractional integrals. 
The order of fractional integral is equal to the fractal mass dimension of the medium. 

A more consistent approach to describe the fractal media is connected with the mathe- 
matical definition of the integrals on fractals. In Ref. [SI], it was proven that integrals on 
net of fractals can be approximated by fractional integrals. In Ref. ^T] we proved that frac- 
tional integrals can be considered as integrals over the space with fractional dimension up to 
numerical factor. To prove this proposition, we use the well-known formulas of dimensional 
regularizations [TU] . 

The fractional continuous models of fractal media ^2] can have a wide application. This 
is due in part to the relatively small numbers of parameters that define a random fractal 
medium of great complexity and rich structure. The fractional continuous model allows 
us to describe dynamics of a wide class of fractal media. Fractional integrals can be used to 
derive the fractional generalization of the equations for the fractal media. 

Let us consider the region W in three-dimensional Euclidean space E 3 . The volume of 
the region W is denoted by V(W). If the region If is a ball with the radius R, then the 
volume V(Wa) = (A/3)ttR a . The mass of the region W in the fractal medium is denoted 
by Md(W), where D is a mass dimension of the medium. 

The properties of the fractal media like mass obeys a power law relation, 



where M is the mass of fractal medium, R is a box size (or a sphere radius), and D is a 
mass fractal dimension. The power law relation ()54|) can be naturally derived by using the 
fractional integral. Let us consider the fractional generalization of the equation 



Let us define the fractional integral in Euclidean space E 3 in the Riesz form (T^j by the 
equation 



M(R) = kR 



,D 



(D<3) 



(54) 




(55) 




(56) 



where dVo = C3(D,r,r )d 3 r, and 



2 3 " D r(3/2) 



3 



c 3 (D,r, r ) 



r — 




r — 
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The point ro € W is the initial point of the fractional integral. We will use the initial points 
in the integrals set to zero (ro = 0). The numerical factor in Eq. (JHH)) has this form in 
order to derive the usual integral in the limit D —>■ (3 — 0). Note that the usual numerical 
factor %\D) = T(l/2)/2 D ir 3 / 2 T(D/2), which is used in Ref. [13], leads to 7^(3 - 0) = 
r(l/2)/2 3 7r 3 / 2 r(3/2) in the limit D -> (3 - 0). 

Using notations (|56|) . we can rewrite Eq. (J55|) in the form M 3 (W) = (I 3 p)(r ). Therefore 
the fractional generalization of this equation can be defined in the form 

M D (W) = (I D p)(r ) = 23 r ^ 2) J w p(r)\r-r \ D -Wr. (57) 

If we consider the homogeneous fractal media (p(r) = p — const) and the ball region 
W, then we have 

2 3 - D r(3/2) 



T(D/2) Jw 

where R = r — r . Using the spherical coordinates, we get 

Md(w) = / R ^ dR = i^m^, 



T(D/2) '"Jw DT(D/2) 

where R = |R|. As the result, we have Eq. JHU), where k = [2 5 - D 7TT(3/2)p ]/[Dr(D/2)} 
Therefore the fractal media with noninteger mass dimension D can be described by fractional 
integral of order D. 
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